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Abstract. We performed a measurement of differential and integral jet shapes in proton- 
carbon, proton-tungsten and proton-aluminium collisions at 920 GeV/c proton momen- 
tum with the HERA-B detector at HERA for the jet transverse energies in the range 

4 < ET(jet) < 12 GeV. Jets were identified using the fcT-clustering algorithm. The mea- 
surements were performed for the hardest jet in the event, directed towards the opposite side 
with respect to the trigger direction. Jets become narrower with increasing transverse energy 
and measured distributions agree well with predictions of the PYTHIA 6.2 model. We do not 
observe any significant difference in the jet shape for the carbon and the aluminium targets. 
Nevertheless, the transverse energy flow at small and large radii for the tungsten sample is 
slightly less than for light nuclei. This observation indicates some influence of the nuclear 
environment on the formation of jets in heavy nuclei, especially at lower transverse energies, 

5 < E T (jet) < 6 GeV. 

PACS. 13.85.-t Hadron-induced high- and super-high-energy interactions - 13. 87. -a Jets in 
large-Q 2 scattering - 13.87.Fh Fragmentation into hadrons 



1 Introduction 

The parton — * hadron transition is one of the most 
interesting outstanding questions in QCD and we 
should look inside jets to better understand this 
process. 

In the Standard Model our general understand- 
ing of high-energy collisions of hadrons suggests that 
jets arise when short-distance, large-momentum-tran- 
sfer interactions generate partons (quarks and glu- 
ons) that are widely separated in momentum space 
just after the hard collision. In a fashion that is not 
yet quantitatively understood in detail these con- 
figurations are thought to evolve into hadronic fi- 
nal states exhibiting collimated sprays of hadrons, 
which are called jets. Thus jets can be regarded as 
a universal signal of parton dynamics at short dis- 
tances. 

Hadron collisions are a perfect place to perform 
such studies, because they are a high rate source 
of jets over a very wide range of QCD scales in 
the same experiment. High-py jets dominate the 
event structure in hadronic collisions at high center- 
of-mass (CM) energies y/s when sufficiently large 
transverse energy Et is required. This was demon- 
strated by various experiments at y/s = 63 GeV 
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HI, at y/s = 540 GeV [2], at y/s = 1800 GeV [3] 
at hadronic colliders. 

At large transverse energies (Et ~ 100 GeV) 
the corrections to jet shape from parton fragmen- 
tation are usually considered to be small in com- 
parison to ones due to the parton cascade initiated 
by the high-E^ scattered parton. But at the more 
moderate transverse energies ~ of a few GeV, this 
conclusion is not valid. Thus the investigation of the 
jet structure at moderate energies in hadronic colli- 
sions can give important information about parton 
fragmentation processes in hadronic interactions. 

Measurements of jets in fixed-target experiments 
not only widen the energy range of jet studies but, 
in addition, enable extension of these studies into a 
new realm of colliding particles: meson- nucleon and 
hadron-nucleus interactions. 

At lower y/s, because of contributions from mech- 
anisms such as initial- and final-state parton ra- 
diation and multiple scattering of quarks and glu- 
ons, the event structure rarely exhibits dijet topol- 
ogy and the jet signal is rather difficult to extract 
experimentally [4j[5]. Nevertheless, the presence of 
jet structure has also been demonstrated in fixed 
target experiments at lower energies specifically at 
800 GeV in pA-collisions 6 , and at 500 GeV in both 
pBe and 7ri?e-collisions 0- Also a jet signal consis- 
tent with QCD predictions and with extrapolations 
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of jet cross sections from higher CM energies was 
found by [HIS]. 

The methods employed to extract jets at mod- 
erate CM energies depend heavily on Monte Carlo 
models. Therefore, the absolute cross sections for jet 
production determined at these energies are subject 
to large systematic errors. However, the relative de- 
pendence of jet production and properties from dif- 
ferent nuclear targets should be less sensitive to the 
assumed jet size and background. 

There is a large variety of jet-shape variables, 
they are very informative and enter all hard pro- 
cesses with features which are expected to be uni- 
versal (QCD factorization). In the present analy- 
sis we concentrate on studies of the internal jet 
structure, measuring the differential and integral jet 
shape. 



2 Jet definition 

A jet is qualitatively defined as a collimated spray 
of high-energy hadrons. However, for the purpose of 
performing accurate quantitative studies, one needs 
a precise definition of a jet. Essentially, one has to 
specify how low-energy particles are assigned to jets, 
in order to have infrared-finite cross sections. 

The standard "Snowmass convention" on jet 3- 
momcntum definition 10J is: 

Er(jet) = J^i PT,t , 

V(jet) = Ei Pt,i Vi /E T {jet) , (1) 
<t>(jet) = J2i Pr,i <f>i I E T (jet) , 

with Exijet), r](jet), 4>(jet) being the jet trans- 
verse energy, pseudo-rapidity and azimuthal angle, 
respectively, and pr,i, Vi? 4>i are transverse energy, 
pseudo-rapidity and azimuthal angle of particles, 
forming the jet. Usage of the above observables en- 
sures the invariance of the jet momentum definition 
with respect to longitudinal Lorentz boosts. 



multi-jet events as well problems when being com- 
pared with NLO and NNLO QCD calculations. 

On the other hand the kx algorithm [TTUT^] is 
expected to be the best for hadronic collisions be- 
cause it is based on the QCD picture of jet develop- 
ment. This algorithm is invariant under longitudinal 
Lorentz boosts and is infrared and collinear safe. It 
has been shown [13] that the inclusive kx cluster 
algorithm provides, at present, the best jet finding 
algorithm from the theoretical point of view, since 
the problem of overlapping jets, which affects, e.g., 
the iterative cone algorithm ,14], is avoided. 

The measure of the "closeness" of two parti- 
cles/protojets in the kr algorithm is: 



dij = min (p TA ,p T 



mf + 



iY 



(2) 



We use here the original expression for dij [TT] 
without the additional parameter Rq suggested in 
[15] . where the measure @ is replaced by dij/R^. 
Because the parameter Rq is to be around unity, it 
does not make sense to introduce Rq in our case. 
At small angles between two particles the measure, 
dij, is, approximately, the squared relative trans- 
verse momentum of one of the particles with respect 
to the other particle. 

The kr algorithm looks as following [TT]: 

1. For each protojet, define di = E\ i and for each 
pair of protojets define dij from ^ 

2. Find the smallest of all the di and dij and label 

it d m i n 

3. If d m in is a dij, merge protojets i and j into a 
new protojet k according to JTJ 

4. If d m i n is a di, the corresponding protojet i is 
"not mergeable" . Remove it from the list of pro- 
tojets and add it to the list of jets. 

Repeating the above four steps till there are 
"mergeable" protojets one subdivides the event into 
a number of groups of particles where each particle 
is assigned to one and only one group/jet. 



Beside the above general definition of the jet's 
characteristics one needs to specify an algorithm 
- how to assign a particle to one of the jets in 
the event. The most usable at present are two al- 
gorithms - "cone" and "fey" (or "Durham") algo- 
rithms. 

The simplest cone algorithm (see, e.g., [TU]) starts 
from a "jet initiator" - a particle with transverse 
momentum above a predefined threshold, e.g., px > 
1 GeV, adjoining to the initiator particles within a 
cone of radius r < R c . Here r — \J Ar/ 2 + A(f> 2 and 
R c ~ 0.7. There are modifications of this algo- 
rithm, using an iterative procedure to define the 
current jet axis. The cone algorithm is suitable for 
e + e~ collisions, where the remnants of the initial 
particles are absent. For hadronic interactions and 
deep-inelastic lepton-hadron scattering, the cone al- 
gorithm meets problems with particles separation in 



3 Jet shape observables 

In the present work we measure the differential p(r) 
and integral tp'(r) jet shapes characterizing how 
widely a jet's energy is spread in the (rj, (j)) plane. Jet 
shape is one of the most popular characteristics of 
jet structure in hadron collisions. Jet shape \P(r) is 
defined as the transverse energy flow within the cone 
of radius r around the jet axis in (rj, (f>) plane nor- 
malized to unity and averaged over all jets in the 
(sub)sample. The differential jet shape p(r) dr is 
the derivative of W(r) over r and is the transverse 
energy flow through the annulus of width dr with 
radius r around the jet axis. 

These observables are, in general, collinear safe 
(for two particles with parallel momenta Pi , Pj one 
mav replace these two particles with a single one 
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with a momentum equal to pi + pj) and infra-red 
safe (for pi <C Pj one may neglect pi) [16] . This im- 
plies that one may assume parton-flow ~ hadron- 
flow. Particles include both charged particles mea- 
sured by the tracker detector and clusters in the 
electromagnetic calorimeter. 

More quantitatively, the observables tf'(r) and 
p{r) are defined as follows: 

*(r)=fodr>p(r>), 

P {r) =2nrJ Q E ^dp T ^ ( 3 ) 

d 2 n(r,p T ) 
dr dpT 5 

where r is the distance to the jet axis in (77, <p) plane, 
Pt is the transverse momentum of a particle, and 

" ' drdp^ * s ^ e P ar ticle number density over r and 
PT- 
In the present work we study also the depen- 
dence of jet structure on atomic number of a target 
nucleus. In principle, there may exist an uncertainty 
with choice of the center of mass reference frame — 
either (pp), or (pn) or (pA). To avoid this problem 
we perform our analysis in the laboratory reference 
frame. 

Due to the above uncertainty with the choice of 
reference frame we do not present the more tradi- 
tional longitudinal and transversal (with the respect 
of the jet axis) distributions of particles ("fragmen- 
tation function"). These observables are not invari- 
ant under longitudinal Lorentz boosts and usually 
are given in the beam-target center of mass system. 



4 HERA-B detector and data sample 
4.1 Detector 

HERA-B is a fixed target experiment operated at 
the 920 GeV proton storage ring of HERA at DESY 

A plan view of the HERA-B spectrometer is 
shown in Fig. [T] The spectrometer dipole magnet 
provides a field integral of 2.13 T-m, with the main 
component perpendicular to the x-z plane. The ap- 
paratus (including particle identification counters) 
has a forward acceptance of 15-220 mrad in the 
bending plane and 15-160 mrad in the non-bending 
plane. The experiment uses a multi-wire fixed target 
which operates in the halo of the proton beam dur- 
ing HERA e-p collider operation. Up to eight dif- 
ferent targets can be operated simultaneously, with 
their positions being adjusted dynamically in order 
to maintain a constant interaction rate between 1 
and 40 MHz. 

The tracking system consists of a Vertex De- 
tector System (VDS) [19] and a main tracker sys- 
tem. The VDS features 64 double-sided silicon mi- 
crostrip detectors arranged in eight stations along, 
and four Quadrants around, the Droton beam. The 



silicon strips have a readout pitch of approximately 
50 fim. Particle identification was performed by a 
Ring Imaging Chercnkov detector (RICH) [22], an 
electromagnetic calorimeter (ECAL) [23] and a muon 
detector (MUON) [M]. 

The main tracking system is separated into an 
Inner Tracker (ITR) [20] close to the proton beam- 
pipe and an Outer Tracker (OTR) [5T] farther out. 
The tracker system covers pseudorapidities within 
a range of approximately 2 < r/ < 4.8 in the 
laboratory system. Magnetic analysis for the OTR 
ensures the relative momentum resolution around 
Ap/p w 5 • 10" 5 p © 1.6 • 10~ 2 . 

The electromagnetic calorimeter ECAL is based 
on "shashlyk" sampling calorimeter technology, con- 
sisting of scintillator layers sandwiched between metal 
absorbers. In the radially innermost section tung- 
sten was used as an absorber, and lead everywhere 
else. ECAL covers pseudorapidities within a range 
of approximately 2 < r/ < 5.2 in laboratory sys- 
tem. Calorimeter towers have a depth of about 20 
radiation lengths with granularity approximately 
2.2 x 2.2 cm, 5.58 x 5.58 cm, and 11.15 x 11.15 cm 
for inner, middle and outer sections respectively. 
Energy resolution AE/E is 20.5%/V^ ffi 1.2% for 

inner section and approximately 11%/Vi? © 1.0% 
for middle and outer sections. 

4.2 Data sample 

For the analysis presented here, only data from VDS, 
OTR and the ECAL were used with carbon (C), alu- 
minum (Al), and tungsten (W) wire targets. The 
inner tracker was not used due to its insufficient 
stability. 

We used data collected in so called "high Et" 
runs, i.e., runs in which about half of the events 
in each run were required to satisfy a calorimeter 
pre-trigger. The "high Et" pre-trigger demands a 
transverse energy deposition at least in one of the 
ECAL towers above a predefined threshold, Et > 
Exirnin). The remaining events in these runs were 
collected with a random interaction trigger, in fact, 
minimum bias events. The data were recorded at a 
moderate interaction rate, 1 — 3 MHz, which corre- 
sponds to ss 0.1 — 0.3 interactions per filled bunch 
crossing. Therefore only a small fraction of events 
contain more than one interaction. Tungsten data 
were collected with pre-scaling by a factor 2. 

The results presented here are based on a sam- 
ple of sa 18 million events collected in the period 
from December 2002 to March 2003. There were 
three runs with a carbon target, two runs with a 
tungsten target, and three runs with an aluminium 
target with ECAL pre-trigger either Et > 3 GeV 
or E T > 2 GeV. Note that about 40% of trig- 
gers were produced by ECAL clusters matched with 
tracks/segments in tracking detector. 

For the additional off-line event selection, we re- 
quire that the ECAL cluster with maximal trans- 
verse momentum (called in the following the "trig- 
ger cluster" ) has vVi ax > 3.0 GeV. The choice of the 
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Fig. 1. Schematic top view of the HERA-B detector 



cluster with maximal pt unifies the selection of data 
and MC samples. This "trigger cluster" does not 
necessarily coincide with the ECAL cluster which 
produced the hardware pre-trigger. But the pre- 
trigger and the hardest ECAL clusters differ in less 
than 0.1% of events. 

Events were required to have a reconstructed 
primary vertex. The track and cluster selection cri- 
teria are the following: 

— Track must start in the VDS, 

— Track must be successfully htted over the whole 
length (OTR and VDS), 

— "Clones" (nearby reconstructed tracks originat- 
ing from the same real physical track) are re- 
moved, 

— EM clusters must not match any track or track 
segment in the OTR, 

— The trigger conditions are tested before the check 
on track-cluster matching. If the trigger cluster 
matches a track or a track segment, the event is 
accepted, but the trigger cluster itself is removed 
from further analysis. 

After applying above selection criteria we have 
for our analysis approximately 1097000 carbon ev- 
ents, 447000 tungsten events and 1002000 aluminium 
events. 

We considered only the jet with the maximal 
transverse energy in the event. Because our trigger 
is based on measuring the electromagnetic clusters 
in the ECAL, it can distort properties of the selected 
jets, preferring the jets enriched by gammas/7r°'s. 
To avoid this problem, we accepted only "away- 
side" jets, i.e., jets directed in the opposite direction 
with respect to the trigger cluster in the plane per- 
pendicular to the beam axis, demanding \<j){jet) — 
cj)(trig)\ > 90°. The Monte Carlo simulation shows 
that such a selection gives minimal distortion of the 
jet transverse energy and its direction in comparison 
to the parent parton. 

To minimize effects of the restricted acceptance 
we accept jets with axes in the narrow pseudo-rapi- 
ditv ranee 3.4 < n(iet) < 3.6. 



5 Monte Carlo simulation 

The simulation of the physical processes and the de- 
tector response in modern experiments is one of the 
crucial parts of data analysis. For the simulation of 
jet events in HERA-B, the main problem is that 
existing generators describing the hadron(nucleus)- 
nucleus collisions, are intended for the study of mini- 
mum bias physics and soft processes. Therefore their 
use for simulation of hard hadron- nucleus collisions 
is extremely inefficient (if it is possible at all) . Due 
to these problems a procedure for the description 
of hard parton scattering, taking into account soft 
interactions in the nucleus was developed. 

The bases of the developed MC generator are 
theoretical expectations as well as experimental evi- 
dence [6j[7] that at moderate energies ^/s « 40 GeV, 
the hard-scattering mechanism already dominates 
pp collisions, whereas soft scattering is still the 
dominant mechanism for the proton collisions with 
heavy nuclei. This interpretation is consistent with 
the decrease of nuclear effects with increasing jet 
Et- All available data on the hadroproduction of 
jets are consistent with a rather modest nuclear en- 
hancement (the parameter a in A-dependence is 
close to unity, a < 1.10), diminishing with increas- 
ing jet E T . 

This fact allows us to discriminate hard and soft 
processes and to exploit the existing MC models for 
hard parton scattering to simulate jet production in 
proton- nucleus collisions. 

We use PYTHIA 6.2 [25] to produce hadronic 
systems in parton-parton scattering, adding to the 
hard system soft particles produced in the collision 
of the proton remnant with the nucleus. Note that in 
PYTHIA, the partonic processes are simulated us- 
ing LO matrix elements with the inclusion of initial- 
and final-state parton showers. 

The program FRITIOF 7.02 [26j from CERN li- 
brary (including ARIADNE 4.02 as a part of FRITI- 
OF package) has been adapted to double precision 
to make it compatible with PYTHIA 6.2. Because 
the FRITIOF package is heavily based on PYTHIA, 
all parameters for soft fsemi-hard") Drocesses in FRI- 
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TIOF are the default parameters of PYTHIA 6.2 in 
the adapted version of FRITIOF. 

In such a combined PYTHIA 6.2/FPJTIOF 7.02 
package the event simulation procedure consists of 
the following steps: 

— simulation of the hard sub-process (es) (parton- 
parton scattering) with PYTHIA in pN collision 
(a nucleon N is being chosen randomly accord- 
ing to the nucleus content), 

— the independent fragmentation of scattered par- 
tons is performed, 

— resulting stable particles (hard system) are saved, 

— simulation of the soft collision is performed by 
the FRITIOF generator as a pA interaction with 
energy remaining after the hard scattering, 

— the transverse momenta of the whole system are 
balanced, preserving the 4-momentum of the hard 
system, 

— the final fragmentation and decays are performed. 

The main PYTHIA inputs used for hard scat- 
tering are the following: 

— Proton structure function "CTEQ2L (best LO 
fit)" from PDFLIB, 

— Intrinsic kr = 1 GeV of partons (default for 
PYTHIA 6.2), 

— Independent fragmentation of final state partons 
from the hard process, 

— String fragmentation in soft interactions with 
the nucleus (in FRITIOF), 

— Lund symmetric fragmentation model both for 
independent and string fragmentation, 

— QCD scale for parton-parton scattering as is the 
default for PYTHIA 6.2. 

We tested that the default set of PYTHIA inputs 
gives a quite satisfactory description of our data as 
well as a good reproduction of the prompt photon 
Pt distribution measured by experiment E-706 [27j 
in pBe collisions at ^/s = 38.8 GeV. 

Note here that because we restricted our studies 
to the narrow range Ar\ = ±0.1 the variation of the 
PDF does not play any role. We also checked that 
the choice of the QCD scale for parton scattering is 
not significant in our range of transverse momenta 
and pseudo-rapidities. 

The main effect comes from the intrinsic trans- 
verse momentum, kx, of partons. In paper |27j the 
authors suggest the value kr = 1.3 GeV to describe 
their results on prompt photon production. But we 
did not find a significant difference between the val- 
ues k T = 1.3 GeV and the default PYTHIA 6.2 
value, kr — 1.0 GeV for our conditions and use the 
default PYTHIA 6.2 value. 

To increase the simulation efficiency the parton 
sub-processes were simulated with the cut 
PT^hard) > 3 GeV, where the px(hard) is the trans- 
verse momentum of the outgoing 7/parton. Such a 
cut increases the simulation efficiency by a few or- 

1 We also tested the simulation of irA collisions for 
the remnant. The difference is insignificant. 



ders of magnitude in comparison to the default cut 
in PYTHIA for hard scattering pT(hard) > 1 GeV. 

The simulation of the detector response was car- 
ried out using the GEANT 3.21 package with subse- 
quent standard HERA-B reconstruction of the sim- 
ulated events. 

We applied the same selection criteria for tracks 
and EM clusters as for real data. In the MC sam- 
ple we also removed from consideration (as in data) 
the so called "hot modules" in ECAL. "Hot mod- 
ules" cover a negligibly small fraction of the total 
ECAL acceptance and, in fact, do not affect the pre- 
sented results. In total, after selection (without trig- 
ger cut) , the number of MC events is approximately 
equal to experimental statistics — 1499000 events 
for the carbon target, 571400 events for the tung- 
sten target and 1016600 events for the aluminium 
target. 



6 Results 

The measured, raw differential jet shape distribu- 
tion p(r) was created by storing the entries in the 
bin Ar with weight pr/ Ex(jet) and further divid- 
ing this distribution by the bin width Ar. 

To correct the measured jet shape distributions 
for acceptance and reconstruction distortions (sys- 
tematic corrections) we used bin-by-bin corrections. 
Corrections for p(r) were calculated in each bin of 
jet transverse energy Et separately and can be 
written as: 

p(r;E T ) =R(r;E T ) x p m (r;E T ) 
R(r;E T )=p g (r;E T )/p r (r;E T ) W 

where, p and p m are corrected and measured distri- 
butions and p r and p g are simulated distributions 
at the reconstructed and generator levels, respec- 
tively. 

We tested the quality of the track information to 
check for possible systematic biases in the track mo- 
mentum measurements. For this purpose we recon- 
structed known resonances K°, p°, K*°, </>°, and A . 
Results for reconstructed masses are in good agree- 
ment with the known values of the Particle Data 
Group [28) . Thus we don't expect any significant 
systematic bias in measured track momenta. 

We checked, that after applying the trigger cut 
p-max > 3 G e v, the simulated distributions for jets 
and particles agree well with data except for the 
multiplicity distributions for low-py particles, where 
simulated events have lower average multiplicity. 
This discrepancy is a known problem of the FRITIOF 
generator. But this discrepancy does not distort jet 
shape distribution in the selected jets and does not 
play a significant role in bin-by-bin correction pro- 
cedure. 

The trigger cut, p™ ax > 3 GeV, leaves about 
0.5% of the total MC sample. Using the simulated 
sample, we checked that the event triggering does 
not distort the iet shaoe for "awav-side" iets. We 
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checked also that triggered and non-triggered MC 
samples both give good description of the data jet 
shape within respective statistical errors. This fact 
allows us to use the non-triggered MC sample to 
find systematic corrections R(r;Er), thus reducing 
statistical systematic errors by more than an order 
of magnitude. 

We performed measurements of the jet shape in 
the range of jet transverse energies Et > 5 GeV, 
to minimize the influence of the cut pxihard) > 
3 GeV which we applied when simulating parton- 
parton scattering. Due to the presence of the Gaus- 
sian-distributed intrinsic transverse momentum with 
width kr = 1 GeV, the missing range of parton 
transverse momenta, pt < 3 GeV, gives a negligi- 
ble contribution to the jet yield for Et > 5 GeV. In 
this case one can expect that MC based systematic 
corrections give reliable results. 

We could select jets up to transverse energy Et < 
14 GeV, however, for Et > 10 GeV, the statistics 
is small due to the steeply falling jet Et spectrum. 
Therefore, we divide Er(jet) range 5—10 GeV into 
five equal bins and chose one bin for Et > 10 GeV. 
We apply bin-by-bin corrections (j4|) to measured 
differential jet shape p(r) in each ET^jet) bin. 

Our results for differential and integral jet shapes 
are presented in Tables [HE] for the carbon target, in 
Tables [3] [4] for the aluminium target and in Tables 
El [6] for the tungsten target. We take as a system- 
atic error the uncertainty on the correction factors 
defined by the finite Monte Carlo statistics. Our 
studies show that MC uncertainties dominate the 
other possible sources of systematic errors. Errors 
given in the tables are the quadratic sum of sta- 
tistical and systematic errors. Both types of errors 
are, approximately, equal each to other. The cor- 
rected ^{t^Et) distribution was obtained by in- 
tegration of the corrected p(t;Et) distribution in 
each Er{jeb) bin according to: 



than jets for light nuclei. However the position of 
maximum energy flow is the same for all nuclear 
targets, r max ps 0.3. 

A jet is not a point-like object and has finite 
transverse size, its radius. In the cone algorithm, 
the jet radius is fixed and equal to the size of the 
cone chosen for jet selection, R cone = 0.7 — 1.0. 
In case of the fc^ algorithm, the jet radius is not 
fixed and extends to values r > 1. In the case of 
restricted acceptance over pseudo-rapidity (in the <fi- 
direction the acceptance is not restricted) one can 
encounter problems with the selection efficiency of 
wide jets, when the jet radius exceeds the detector 
acceptance. We performed measurements for jets 
with axes within the pseudo-rapidities r/(jet) — 3.5± 
0.1. A realistic acceptance of HERA-B tracker (prob- 
ability to register a track is greater than 10%) covers 
the pseudo-rapidity range 2.2 < rj < 4.8 (calorime- 
ter has larger acceptance). Therefore, some fraction 
of selected jets must have particles outside the ac- 
ceptance of our detector. However, the integral jet 
shape i^(r) is about unity already for r f» 1, i.e., al- 
most all of the jet energy flow is concentrated within 
this cone and the whole jet is within the detector 
acceptance (2.2 < r) < 4.8). 

Fig|4]demonstrates the correction factor R(r; Et) 
from ([4]) normalized to unity for the carbon tar- 
get in the bin E T {jet) = 5.0 - 6.0 GeV. The hor- 
izontal line is the mean value for this distribution. 
The variation of the correction factor does not ex- 
ceed ~ 25% with respect to its mean value at radii 
r < 1.5 with good statistical errors. Thus we can 
expect that our systematic corrections are reliable 
within the effective interval of jet radii. Note here 
that correction factors for different nuclear targets 
agree well within their errors. Therefore they cancel 
out in the ratios of p(r) for different nuclei, pre- 
sented in Fig 03 



8=1 

where i denotes the bin in a histogram. 

Results of measurements of p[r) and W{r) are 
presented in Tables [THE] 

Fig[2] shows the comparison of the differential 
jet shape p{r) for carbon target for six bins of jet 
transverse energy (closed circles) and PYTHIA/- 
FRITIOF predictions (open circles) . The agreement 
between data and model predictions is quite good. 

Fig[3] presents the ratios of the p(r) measured 
with aluminium (closed circles) and tungsten (open 
circles) targets to carbon results. The results for 
aluminium and carbon targets are nearly the same 
(ratio is approximately equal to unity) . Statistically 
significant deviations from unity exist for the tung- 
sten sample at r w in Et range 5 < Et < 
10 GeV and for large radius, r > 1 in the bin 
5 < Et < 6 GeV. Thus, on average, jets produced 
on tungsten nuclei have sliehtlv different structure 
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Fig. 4. The dependence of the systematic correction 
factor R(r; Et) on the distance to the jet axis for carbon 
target in Erijet) — 5.0 — 6.0 GeV bin. The correction 
factor is normalized to unity. 
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Table 1. Differential and integral jet shapes for carbon for 5 < Et < 8 GeV. 





E T = 5-6 GeV 












Et z 


= 6- 


- 7 GeV 












Et - 


= 7- 


- J- 


GeV 


r 










p 












p 




















p 












u 


uo 





550 


± 


o 


040 


055 ± 004 





550 


± 





050 





055 


± 





005 





420 


± 





050 





042 


± 





005 


n 
u 


1 K 
10 


1 


250 


± 


o 


050 


1 SO + 006 


1 


210 


± 





060 


o 


176 


± 


o 


008 


1 


050 


± 





080 


o 


147 


± 


o 


009 


n 
u 


oc: 
ZO 


1 


510 


± 


o 


050 


331 ± 008 


1 


390 


± 


o 


060 


o 


320 


± 


o 


010 


1 


500 


± 


o 


100 


o 


295 


± 


o 


014 


n 
u 


OO 


1 


440 


± 


o 


040 


476 ± 009 


1 


410 


± 





050 





460 


± 


o 


010 


1 


430 


± 





080 





438 


± 


o 


016 


n 
u 




1 


370 


± 


o 


040 


613 ± 009 


1 


410 


± 


o 


050 





600 


± 


o 


010 


1 


420 


± 





070 


o 


580 


± 


o 


020 


n 
u 


00 


1 


200 


± 


o 


030 


730 ± 010 


1 


120 


± 





040 





710 


± 


o 


010 


1 


120 


± 





060 





690 


± 


o 


020 





65 


o 


880 


± 


o 


020 


820 ± 010 





890 


± 





030 





800 


± 


o 


010 





900 


± 





040 





780 


± 


o 


020 





75 





630 


± 





020 


0.880 ± 0.010 


o 


680 


± 





020 





866 


± 





014 


o 


750 


± 





040 





860 


± 





020 





85 





460 


± 





010 


0.930 ± 0.010 





490 


± 





020 





915 


± 





014 





470 


± 





020 





900 


± 





020 





95 





287 


± 





009 


0.960 ± 0.010 





358 


± 





014 





951 


± 





014 





380 


± 





020 





940 


± 





020 


1 


05 





177 


± 





006 


0.980 ± 0.010 





220 


± 





010 





973 


± 





014 





242 


± 





015 





970 


± 





020 


1 


15 





112 


± 





005 


0.990 ± 0.010 





122 


± 





006 





985 


± 





014 





150 


± 





010 





980 


± 





020 


1 


25 





057 


± 





003 


0.990 ± 0.010 





076 


± 





004 





992 


± 





014 





097 


± 





008 





990 


± 





020 


1 


35 





032 


± 





002 


1.000 ± 0.010 





039 


± 





003 





996 


± 





014 





049 


± 





005 


1 


000 


± 





020 


1 


45 





017 


± 





001 







019 


± 





002 





998 


± 





014 





024 


± 





003 












1 


55 





008 


± 





001 







011 


± 





001 





999 


± 





014 





011 


± 





002 












1 


65 





003 


± 





001 







005 


± 





001 


1 


000 


± 





014 





006 


± 





001 












1 


75 





001 


± 





001 







001 


± 





001 















003 


± 





001 












1 


85 

















001 


± 





001 















001 


± 





001 













Table 2. Differential and integral jet shapes for carbon for Et > 8 GeV. 















Et 1 


= 8- 


-9 GeV 


E T = 9 - 10 GeV 


E T > 10 GeV 


r 










p 




















p 












p 















05 





280 


± 





040 





028 


± 





004 





280 


± 





050 


0.027 ± 0.005 





150 


± 





060 





015 


± 





006 





15 





900 


± 





100 





121 


± 





010 





700 


± 





100 


0.100 ± 0.010 


1 


000 


± 





300 





120 


± 





030 





25 


1 


350 


± 





140 





256 


± 





020 


1 


140 


± 





160 


0.220 ± 0.020 


1 


600 


± 





300 





280 


± 





040 





35 


1 


300 


± 





100 





383 


± 





020 


1 


600 


± 





200 


0.370 ± 0.030 


1 


400 


± 





300 





420 


± 





050 





45 


1 


400 


± 





100 





529 


± 





020 


1 


600 


± 





200 


0.530 ± 0.040 


1 


400 


± 





200 





560 


± 





050 





55 


1 


220 


± 





090 





651 


± 





030 


1 


200 


± 





100 


0.650 ± 0.040 


1 


130 


± 





150 





670 


± 





060 





65 


1 


040 


± 





070 





755 


± 





030 





960 


± 





090 


0.750 ± 0.040 





670 


± 





090 





740 


± 





060 





75 





790 


± 





060 





835 


± 





030 





650 


± 





060 


0.810 ± 0.040 





640 


± 





090 





800 


± 





060 





85 





550 


± 





040 





889 


± 





030 





620 


± 





060 


0.880 ± 0.040 





560 


± 





080 





860 


± 





060 





95 





430 


± 





030 





933 


± 





030 





380 


± 





040 


0.910 ± 0.040 





470 


± 





080 





910 


± 





060 


1 


05 





260 


± 





020 





959 


± 





030 





340 


± 





040 


0.950 ± 0.040 





360 


± 





060 





940 


± 





060 


1 


15 





200 


± 





020 





979 


± 





030 





220 


± 





030 


0.970 ± 0.040 





130 


± 





030 





960 


± 





060 


1 


25 





100 


± 





010 





989 


± 





030 





130 


± 





020 


0.980 ± 0.040 





190 


± 





040 





970 


± 





060 


1 


35 





050 


± 





007 





994 


± 





030 





084 


± 





015 


0.990 ± 0.040 





130 


± 





030 





990 


± 





060 


1 


45 





023 


± 





004 





996 


± 





030 





050 


± 





010 


1.000 ± 0.040 





050 


± 





020 





990 


± 





060 


1 


55 





015 


± 





003 





998 


± 





030 





028 


± 





009 







040 


± 





010 


1 


000 


± 





060 


1 


65 





015 


± 





006 





999 


± 





030 





015 


± 





006 







040 


± 





020 












1 


75 





002 


± 





001 


1 


000 


± 





030 





005 


± 





002 







007 


± 





004 












1 


85 





002 


± 





002 















004 


± 





003 







008 


± 





005 












1 


95 





































003 


± 





003 













7 Conclusion 



We performed measurements of the differential p(r) 
and integral ^(r) jet shapes ([3]) in the range of 
jet transverse energies 5 < Exijet) < 14 GeV 
in proton-nucleus collisions at a proton momentum 
920 GeV/c. For jet selection, we used the longitudi- 
nally invariant kx algorithm, defining jets according 
to standard Snowmass convention. We performed 
these measurements for three target nuclei — car- 
bon, aluminium and tungsten. For systematic cor- 
rections to the differential jet shape we used bin-by- 
bin statistical corrections, based on a Monte Carlo 
simulation of ict production, in each bin of iet trans- 



verse energy. Integral jet shapes have been obtained 
by integration of the differential jet shapes. 

We find good agreement between data and the 
predictions of the PYTHIA 6.2/FRITIOF 7.02 model 
with intrinsic transverse momentum of partons kr = 
1 GeV and independent fragmentation of partons 
according to the symmetric Lund scheme. The had- 
ronization scheme in PYTHIA 6.2 at studied jet 
transverse energies does not include showering of 
the scattered parton. It instead produces the final 
state hadrons, mostly vector mesons, directly from 
the scattered parton. The agreement between our 
measurements and simulation shows that in the con- 
sidered transverse enerev ranee E-r(iet) < 14 GeV. 
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Table 3. Differential and integral jet shapes for aluminium 5 < Et < 8 GeV. 





E T = 5-6 GeV 








E T = 6 - 


-7 GeV 








Et = 7 - 


-8 GeV 


r 






p 


<f 






P 












P 






0.05 


0.660 


± 


0.050 


066 ± 005 


0.580 


± 


0.060 


0.058 


± 


0.006 


0.430 


± 


0.050 


0.043 


± 0.005 


0.15 


1.230 


± 


0.050 


189 + 007 


1.220 


± 


0.070 


0.179 


± 


0.009 


1.160 


± 


0.100 


0.160 


± 0.010 


0.25 


1.480 


± 


0.050 


337 ± 009 


1.510 


± 


0.080 


0.330 


± 


0.010 


1.400 


± 


0.100 


0.303 


± 0.016 


0.35 


1.450 


± 


0.040 


480 ± 010 


1.420 


± 


0.060 


0.470 


± 


0.010 


1.480 


± 


0.090 


0.450 


± 0.020 


0.45 


1.380 


± 


0.040 


620 ± 010 


1.340 


± 


0.050 


0.610 


± 


0.010 


1.250 


± 


0.070 


0.580 


± 0.020 


0.55 


1.120 


± 


0.030 


730 ± 010 


1.150 


± 


0.040 


0.722 


± 


0.015 


1.060 


± 


0.060 


0.680 


± 0.020 


0.65 


0.870 


± 


0.020 


820 ± 010 


0.850 


± 


0.030 


0.806 


± 


0.015 


0.950 


± 


0.050 


0.780 


± 0.020 


0.75 


0.660 


± 


0.020 


880 + 010 


0.670 


± 


0.020 


0.874 


± 


0.016 


0.750 


± 


0.040 


0.850 


± 0.020 


0.85 


0.450 


± 


0.014 


0.930 ± 0.010 


0.480 


± 


0.020 


0.922 


± 


0.016 


0.560 


± 


0.030 


0.910 


± 0.020 


0.95 


0.290 


± 


0.010 


0.960 ± 0.010 


0.310 


± 


0.010 


0.953 


± 


0.016 


0.350 


± 


0.020 


0.940 


± 0.020 


1.05 


0.194 


± 


0.008 


0.980 ± 0.010 


0.200 


± 


0.010 


0.973 


± 


0.016 


0.247 


± 


0.016 


0.970 


± 0.020 


1.15 


0.111 


± 


0.005 


0.990 ± 0.010 


0.123 


± 


0.007 


0.985 


± 


0.016 


0.150 


± 


0.010 


0.980 


± 0.020 


1.25 


0.063 


± 


0.003 


0.990 ± 0.010 


0.074 


± 


0.005 


0.993 


± 


0.016 


0.086 


± 


0.007 


0.990 


± 0.020 


1.35 


0.030 


± 


0.002 


1.000 ± 0.010 


0.038 


± 


0.003 


0.996 


± 


0.016 


0.042 


± 


0.005 


1.000 


± 0.020 


1.45 


0.016 


± 


0.001 




0.022 


± 


0.002 


0.999 


± 


0.016 


0.024 


± 


0.003 






1.55 


0.008 


± 


0.001 




0.009 


± 


0.001 


0.999 


± 


0.016 


0.009 


± 


0.002 






1.65 










0.004 


± 


0.001 


1.000 


± 


0.016 


0.009 


± 


0.002 






1.75 










0.001 


± 


0.001 








0.004 


± 


0.001 






1.85 






















0.002 


± 


0.001 






1.95 






















0.001 


± 


0.001 







Table 4. Differential and integral jet shapes for aluminium Et > 8 GeV. 











Et — 8 - 


-9 GeV 








E T = 9 - 


10 GeV 


E T > 10 GeV 


r 






P 






<f 






P 












P 








0.05 


0.340 


± 


0.050 


0.034 


± 


0.005 


0.300 


± 


0.060 


0.030 


± 


0.006 


0.170 


± 


0.060 


0.017 


± 


0.006 


0.15 


1.000 


± 


0.100 


0.130 


± 


0.010 


0.600 


± 


0.100 


0.090 


± 


0.010 


0.900 


± 


0.200 


0.110 


± 


0.020 


0.25 


1.440 


± 


0.150 


0.280 


± 


0.020 


1.300 


± 


0.200 


0.220 


± 


0.020 


1.200 


± 


0.200 


0.230 


± 


0.030 


0.35 


1.300 


± 


0.100 


0.400 


± 


0.020 


1.300 


± 


0.200 


0.360 


± 


0.030 


1.600 


± 


0.300 


0.390 


± 


0.040 


0.45 


1.300 


± 


0.100 


0.530 


± 


0.030 


1.300 


± 


0.150 


0.490 


± 


0.030 


1.000 


± 


0.200 


0.490 


± 


0.040 


0.55 


1.300 


± 


0.100 


0.660 


± 


0.030 


1.300 


± 


0.100 


0.610 


± 


0.040 


1.200 


± 


0.200 


0.610 


± 


0.050 


0.65 


1.030 


± 


0.080 


0.760 


± 


0.030 


1.000 


± 


0.100 


0.710 


± 


0.040 


0.900 


± 


0.100 


0.700 


± 


0.050 


0.75 


0.740 


± 


0.050 


0.830 


± 


0.030 


0.900 


± 


0.100 


0.800 


± 


0.040 


0.550 


± 


0.080 


0.750 


± 


0.050 


0.85 


0.500 


± 


0.040 


0.880 


± 


0.030 


0.600 


± 


0.060 


0.860 


± 


0.040 


0.700 


± 


0.200 


0.820 


± 


0.060 


0.95 


0.480 


± 


0.040 


0.930 


± 


0.030 


0.540 


± 


0.060 


0.920 


± 


0.040 


0.600 


± 


0.100 


0.880 


± 


0.060 


1.05 


0.270 


± 


0.020 


0.960 


± 


0.030 


0.330 


± 


0.040 


0.950 


± 


0.040 


0.400 


± 


0.100 


0.920 


± 


0.060 


1.15 


0.152 


± 


0.015 


0.970 


± 


0.030 


0.180 


± 


0.030 


0.970 


± 


0.040 


0.400 


± 


0.100 


0.960 


± 


0.060 


1.25 


0.121 


± 


0.014 


0.990 


± 


0.030 


0.160 


± 


0.030 


0.980 


± 


0.040 


0.180 


± 


0.030 


0.980 


± 


0.060 


1.35 


0.070 


± 


0.010 


0.990 


± 


0.030 


0.082 


± 


0.016 


0.990 


± 


0.040 


0.100 


± 


0.020 


0.990 


± 


0.060 


1.45 


0.037 


± 


0.007 


1.000 


± 


0.030 


0.048 


± 


0.010 


1.000 


± 


0.040 


0.060 


± 


0.020 


1.000 


± 


0.060 


1.55 


0.024 


± 


0.006 








0.025 


± 


0.007 








0.030 


± 


0.010 








1.65 


0.008 


± 


0.002 








0.010 


± 


0.004 








0.008 


± 


0.003 








1.75 


0.003 


± 


0.001 








0.001 


± 


0.001 








0.007 


± 


0.003 








1.85 


0.005 


± 


0.003 
































1.95 


0.001 


± 


0.001 

































the parton cascade is still not developed and the 
hadronization has mainly a non-perturbative na- 
ture, i.e., direct transition parton — > hadrons and 
can be well reproduced by PYTHIA 6.2 with de- 
fault parameters and independent fragmentation of 
outgoing partons. 

From comparison of the differential jet shape 
p(r) for different nuclei, we can conclude that the 
differences in jet properties for carbon and aluminium 
targets are small. Nevertheless, the transverse en- 
ergy flow at small and large radii for the tungsten 
sample is slightly less than for light nuclei. Note 
here that this observation does not depend on our 
MC model due to cancellation of the systematic 



corrections in the ratios of p{r) for different nu- 
clei. This observation indicates some influence of 
the nuclear environment on the formation of jets in 
heavy nuclei, especially at lower transverse energies, 
5 < Exijet) < 6 GeV. Possibly such a "broaden- 
ing" of jets produced on a heavy nucleus occurs due 
to re-scattering of the jet hadrons on nucleons which 
effect must be small for light nuclei. 



We express our gratitude to the DESY laboratory for 
the financial and technical support during this work. 
We are also indebted to the HERA-B Collaboration for 
permission to use the experimental data and for fruit- 
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Table 5. Differential and integral jet shapes for tungsten 5 < Et < 8 GeV. 















Et '- 


= 5- 


- ( 


GeV 


E T = G-7 GeV 












Et - 


= 7- 


- J- 


GeV 


r 










p 




















p 












p 












u 


uo 





490 


± 





050 





049 


± 





005 





390 


± 





040 


039 ± 004 





280 


± 





040 





028 


± 





004 


n 

U 


1 K 
10 


1 


040 


± 


o 


060 


o 


154 


± 





008 


1 


030 


± 


o 


070 


142 + 009 


o 


790 


± 





070 


o 


107 


± 





008 


n 
u 


zo 


1 


390 


± 


o 


060 


o 


290 


± 


o 


010 


1 


380 


± 


o 


080 


280 ± 010 


1 


200 


± 


o 


100 


o 


230 


± 


o 


010 


n 
u 


OO 


1 


580 


± 





060 





450 


± 


o 


010 


1 


460 


± 


o 


070 


430 ± 010 


1 


460 


± 





090 





374 


± 


o 


016 


n 
u 




1 


420 


± 


o 


050 


o 


590 


± 


o 


010 


1 


430 


± 


o 


060 


569 ± 015 


1 


570 


± 





090 


o 


530 


± 


o 


020 


n 
u 


00 


1 


280 


± 





040 





720 


± 


o 


010 


1 


220 


± 


o 


050 


691 ± 016 


1 


320 


± 





070 





660 


± 


o 


020 





65 


1 


000 


± 





040 





820 


± 


o 


010 


1 


080 


± 


o 


040 


800 ± 020 


1 


180 


± 





060 





780 


± 


o 


020 





75 





730 


± 


o 


030 


o 


894 


± 





015 





760 


± 





030 


0.880 ± 0.020 


o 


780 


± 





040 


o 


860 


± 





020 





85 





460 


± 





020 





940 


± 





015 





510 


± 





020 


0.930 ± 0.020 





550 


± 





030 





910 


± 





020 





95 





280 


± 





010 





968 


± 





015 





340 


± 





020 


0.960 ± 0.020 





370 


± 





020 





950 


± 





020 


1 


05 





154 


± 





008 





983 


± 





015 





180 


± 





010 


0.980 ± 0.020 





223 


± 





016 





970 


± 





020 


1 


15 





088 


± 





005 





992 


± 





015 





114 


± 





008 


0.990 ± 0.020 





140 


± 





010 





990 


± 





020 


1 


25 





043 


± 





003 





996 


± 





015 





046 


± 





004 


1.000 ± 0.020 





067 


± 





007 





990 


± 





020 


1 


35 





022 


± 





002 





998 


± 





015 





031 


± 





003 







037 


± 





005 


1 


000 


± 





020 


1 


45 





009 


± 





001 





999 


± 





015 





014 


± 





002 







010 


± 





002 












1 


55 





004 


± 





001 


1 


000 


± 





015 





005 


± 





001 







008 


± 





002 












1 


65 





002 


± 





001 















002 


± 





001 







002 


± 





001 












1 


75 





001 


± 





001 















001 


± 





001 







002 


± 





001 












1 


85 





































001 


± 





001 













Table 6. Differential and integral jet shapes for tungsten Et > 8 GeV. 





E T = 8-9 GeV 












Et 


9- 


10 GeV 


E T > 10 GeV 


r 










p 












p 




















p 















05 





180 


± 
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Fig. 3. Ratios of the p(r) measured with aluminium (closed circles) and tungsten (open circles) targets to carbon 
results in every Er(jet) bin. 



